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Abstract 

The dynamical nucleus-nucleus potentials for fusion reactions 40 Ca+ 40 Ca, 48 Ca+ 208 Pb and 
126 Sn+ 130 Te are studied with the improved quantum molecular dynamics (ImQMD) model to- 
gether with the extended Thomas-Fermi approximation for the kinetic energies of nuclei. The 
obtained fusion barrier for 40 Ca+ 40 Ca is in good agreement with the extracted fusion barrier from 
the measured fusion excitation function, and the depth of the fusion pockets are close to the re- 
sults of time-dependent Hartree-Fock calculations. The energy dependence of fusion barrier is also 
investigated. For heavy fusion system, the fusion pocket becomes shallow and almost disappears 
for symmetric systems and the obtained potential at short distances is higher than the adiabatic 
potential. 
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I. INTRODUCTION 



The synthesis of super-heavy elements (SHEs) has been studied for many years both 
theoretically and experimentally flfl. Up to now the super-heavy nuclei are uniquely 
synthesized through fusion reaction including "cold" fusion reaction with lead and bismuth 
targets [3] and "hot" fusion with actinide targets 4]. The study of the dynamical process 
in fusion reactions especially the nucleus-nucleus potential is of great importance for the 
synthesis of SHEs. Experimentally, the fusion barrier distributions can be directly obtained 
from the measured fusion excitation functions, with which the information of the nucleus- 
nucleus potential around the fusion barrier can be obtained. Fig. 1(a) shows the nucleus- 
nucleus potential calculated with different models for 40 Ca+ 40 Ca. We can see that the 
obtained barrier heights with different models are close to each other and all of them are 
comparable with the extracted mean barrier height, while the calculated nucleus- nucleus 
potentials at short distances are quite different with different models. It is known that the 
adiabatic and diabatic approximations lead to different nucleus-nucleus potentials especially 
at short distances and thus to different fusion paths and different mechanisms of fusion 
reactions. Both approximations are frequently applied to the study of the synthesis of 
super-heavy nuclei [lO] • For understanding the fusion mechanism of a heavy system, it 
is important and necessary to study the nucleus-nucleus potential at short distances, with 
which one could get information on the fusion path and the formation probability of the 
di-nuclear system in reactions leading to super-heavy nuclei. 

For description of heavy-ion fusion reactions, some theoretical models have been devel- 
oped. The fusion coupled channel model is a powerful tool to calculate the fusion excitation 
function and to investigate the influence of nuclear structure effects on the fusion cross 



sections 
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13] . On the other hand, the microscopic dynamics model such as time depen- 
dent Hartree Fock (TDHF) model Q, [l5| and the improved quantum molecular dynamics 
(ImQMD) model Q, Q are widely applied to study the dynamical behavior of fusion pro- 
cess. The ImQMD model is a semi-classical microscopic dynamics model and is successfully 



used for intermediate-energy 



near the Coulomb barrier 



reavy-ion collisions and for heavy-ion collisions at energies 



16Hl9j]. In the ImQMD model the dynamical effects such as the 



dynamical deformation, neck formation, etc, are microscopically and self-consistently taken 
into account. Recently, the model has been applied to the study of the dynamical barrier 
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FIG. 1: (Color online) (a) The nucleus- nucleus potential as a function of the center-to-center 
distance R between two nuclei for the reaction Ca+ 40 Ca. The solid curve denotes the results of 
Skyrme energy density functional together with the extended Thomas-Fermi (ETF2) approximation 



5]. The dash-dotted curve and the crosses denote the Bass potential 



6] and the proximity potential 



7], respectively, (b) The fusion barrier distribution extracted from the measured fusion excitation 
function jg] (see Eq.(9) of Ref.fl]). 



of a heavy 



238 



system 



22 



20], of mass parameters 2l|, and of the strongly damped process of 



23] . In this paper we carefully investigate the kinetic energies of nuclei based 



on the extended Thomas-Fermi approximation with which the dynamical fusion barrier is 
accurately obtained. The paper is organized as follows: In Sec. II, the ImQMD model 
will be briefly introduced. In Sec. Ill, some calculated results on the kinetic energies of 
nuclei and the nucleus-nucleus potential for 40 Ca+ 40 Ca, 48 Ca+ 208 Pb and 126 Sn+ 130 Te will 
be presented. Finally, a summary is given in Sec. IV. 



II. THE IMPROVED QUANTUM DYNAMICS MODEL 



In the ImQMD model, the same as in the original QMD model 24j, each nucleon is rep- 
resented by a coherent state of a Gaussian wave packet. Through a Wigner transformation, 
one can get the one-body phase space distribution function for N-distinguishable particles 
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(see jig, |24j for details). The density distribution function p of a system reads 

\2 

(2^)3/2 C " P[ 2of 

Where, a r represents the spatial spread of the wave packet. The propagation of nucleons 
is governed by Hamiltonian equations of motion under the self-consistently generated mean 
field: 

dH . dH 
r ' = 5p-' P ' = -^- (2) 
Where, and are the centers of i-th wave packet in the coordinate and momentum space, 
respectively. The Hamiltonian H consists of the kinetic energy and the effective interaction 
potential energy: 

H = T + U, (3) 

r = y^. (4) 

i 

The effective interaction potential energy includes the nuclear interaction potential energy 
and the Coulomb interaction potential energy, 

U = U loc + U Couh (5) 

with 



Where Vi oc (r) is the potential energy density w 
interaction and taken as the same as that in Ref. 



U loc = I V loc (r)dr. (6) 
rich is obtained by the effective Skyrme 



V loc = ^ + ^-C^ + f^(Vp) 2 + g T ^ + ^(p 2 - n s {Vpf)5\ (7) 
2 po 7 + 1 p ' 2p Po' 2p 

where 5 = p p n + p p p is the isospin asymmetry. Inserting expression (1) together with (7) into 
(6), we obtain the interaction potential energy 

Pi (3 sr^ f Pi V ! 9o y^y^ f Pij 



Uloc = ^y^ + ^y « ( 8 ) 

2^-po 1 + l^\pJ Z^ 2 -' Po KJ 



i j^i 



where 



Pi = E pa = E jj^m exp [-^^]' (9) 



TABLE I: Parameters set IQ2. 



Parameter 


a(MeV) /3(MeV) 


7 


5o (MeVfm 2 ) g T {MeV) 


rj 


C s (MeV) 


K s (fm 2 ) p ({m 3 ) 


IQ2 


-356 303 


7/6 


7.0 12.5 


2/3 


32.0 


0.08 0.165 



2<x 2 



2a 2 r 



(10) 



and ti = 1 for protons and —1 for neutrons respectively. The parameters set IQ2 [It 
(see Table 1) is adopted in this work. The Coulomb energy is written as the sum of the 



direct and the exc 
approximation 25 



range contribution, and the latter being taken into account in the Slater 



26| 



U, 



Coul 



(11) 



To describe the fermionic nature of the N-body system and to improve the stability of 
an individual nucleus, the phase space occupation constraint method 27J and the system- 



size-dependent wave-packet width a r = 0.09A 1 / 3 + 0.88 fm [18] are adopted. The phase 
space occupation constraint is an effective approach to improve the momentum distribution 
of nuclear system Q- I n this approach, the phase space occupation number of each 
particle is checked at each time step. If the phase space occupation number is larger than 
1 for particle i, i.e. /j > 1, the momentum of the particle i are randomly changed by a 
series of two-body elastic scattering between i and its partner which guarantee that the 
total momentum and total kinetic energy are conserved in the procedures. In the ImQMD 
model, the new sample for the momenta of the particles is constrained by the Pauli-blocking 
probability [27] as in the usual two-body collision process. Actually, the momenta of two 
particles obtained in this way not only influence the motion of particles of the system in 
this step but also the further more steps. It is unknown whether the system will be in the 
most suitable motion path. In this work, we perform one more step further, i.e. we calculate 
the total energy of the system at step t and also the total energy E(t + At) at the next 
time step (t + At) simultaneously. If the value of E(t + At) obviously deviates from that 
of E(t), the two-body elastic scattering procedure is re-executed. The number of times of 
the re-executing process is small (about ~ 4) at each time step for fusion reactions. This 
additional constraint can further improve the stability of an individual nucleus (reducing 



the spurious emission of nucleons), and is helpful for the study of the formation process of 
the compound nuclei which lasts several thousand fm/c or longer. We have checked that 
the total energy of system is well conserved for thousands of fm/c with this new procedure. 

III. RESULTS 

In this section we first study the kinetic energies of a series of nuclei. Then, we calcu- 
late the nucleus-nucleus potential in fusion reactions based on the extended Thomas-Fermi 
approximation. 

A. Kinetic energies of nuclei 

We first study the kinetic energy of a series of nuclei in the ground state from 16 O to 
259 No. Based on the extended Thomas- Fermi (ETF) approximation [28], the kinetic energy 
of a free Fermi gas can be expressed as 

Ek= ^j T{I)dI= ^ + «Li_jVf>: dI+ ... (12) 

with the kinetic energy density r(r) and the coefficient = ^~|(^f~) 2 ^ 3 - With the help of 
the ETF form of the kinetic energy for Fermi gas system in Eq.(12), we express the kinetic 
energy of an individual nucleus in the ImQMD model as 

- co £ P T + jp- £ + ( 13 ) 

with c = 41.2 MeVfm 2 , c\ = 4.8 MeVfm 2 and c 2 = —1.0 MeV for IQ2, which are determined 
by fitting the obtained kinetic energies of a series of nuclei with Eq.(4), see Fig. 2. N is the 
particle number of the system under consideration. The expressions of p$, pij and f s are 
previously given in Eq.(9) and (10), respectively. The cq term of Eq.(13) represents the 
result of the Thomas-Fermi approximation (see the (p) 2//3 term of Eq.(12)). The other terms 
give the corrections from the finite system effect. 

We show the time average of the kinetic energy per particle for a series of nuclei in Fig. 
2(a). Here we take 100 events for each nucleus. The solid and open circles denote the 
results with Eq.(4) and Eq.(13) which is based on the "ETF" approximation, respectively. 
Here, "ETF" means that the form of Eq.(13) is roughly obtained according to the extended 
Thomas- Fermi approximation. The crosses denote the results with the "TF" approximation 

£ fc TF -co$> 2/3 , (14) 
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FIG. 2: (Color online) (a) Time average of the kinetic energy per particle for a series of nuclei. 
The solid circle, the open circle and the crosses denote the results with Eq.(4), (13) and (14), 
respectively, (b) The time evolution of the kinetic energy per particle for 40 Ca, 90 Zr, and 208 Pb. 

in which the correction terms from the finite system effect are not taken into account. From 
the figure one can see that for the light nuclei (A < 50), the calculated kinetic energies 
with the "TF" approximation are much smaller than the values with Eq.(4), while for heavy 
nuclei (A > 150), the results with the "TF" approximation are slightly larger than those 
with Eq.(4). Only for the intermediate nuclei the results are in agreement with each other. 
The calculated kinetic energy with the " ETF" approximation is in good agreement with the 
values with Eq.(4) except for very light nuclei. From Fig. 2, one can see that the "ETF" 
approximation can reasonably well describe the kinetic energy for finite nuclei. The time 
evolution of the kinetic energies per particle for nuclei 40 Ca, 90 Zr, and 208 Pb are shown in 
Fig. 2(b). The kinetic energies of these nuclei can be well described by the Eq.(13) based 
on the "ETF" approximation. 

B. Nucleus- nucleus potential in fusion reactions 

By using the ImQMD model, we can calculate the static and dynamical Coulomb barri- 
ers. In the calculation of the static Coulomb barrier which is based on the frozen density 
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approximation, the initial density distribution of the projectile and target is adopted. In 
the calculation of the dynamical Coulomb barrier, the realistic density distribution of the 
system which changes with time due to the interaction between nucleons is used. In this 
work, we concentrate on the calculation of the dynamical fusion barrier. We study the 
dynamical nucleus-nucleus potential V based on the "ETF" approximation for the kinetic 
energy. According to the energy conservation, we have 

E c .m. = T R + V + E* + T oth , (15) 

where E c , m , is the incident center-of-mass energy, Tr is the relative motion kinetic energy 
of two colliding nuclei, which can be easily obtained in the ImQMD model since the po- 
sition and momentum of each nucleon can be followed at every time step in this model, 
E* is the excitation energy, T ot h is other collective kinetic energy, such as neck vibration. 
When the projectile and target nucleus are well separated (R R\ + R 2 ), the E* and T ot h 



could be neg! 
calculations 



igible which have been checked by the time-dependent Hartree-Fock (TDHF) 



29], the nucleus- nucleus potential is thus expressed as 

Vi = E c . m . - T R . (16) 
Where, R± and R 2 are the charge radii of the projectile and the target nucleus, respectively, 



which are described by an empirical formula i?, = 1.25A 1//3 (1 — 0.2 N . Z ) proposed in 30 ]. 



After the di-nuclear system is formed (R < Ri + R 2 ), the nucleus-nucleus potential may be 



described by a way like the entrance channel potential 31 ] 



V 2 = E tat (R)-E 1 -E 2 , (17) 

where E tot (R) is the energy of the composite system which is strongly dependent on the 
dynamical density distribution of the system obtained with the ImQMD model, Ei and E 2 
are the time average of the energies of the projectile and target nuclei, respectively. Here, 
the values of E\ and E 2 are obtained from the energies of the projectile (like) and target 
(like) nuclei in the region R T < R < R T + 8. R? = Ri + R 2 is the touching point. R is the 
relative distance between the two nuclei, which is a function of time. In the calculation of 
Etot(R), Ei and E 2 , Eq.(13) that is a function of local density is used for the description of 
the intrinsic kinetic energy of the system under consideration. 
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FIG. 3: (Color online) The dynamical nucleus-nucleus potential of 40 Ca+ 40 Ca at different incident 
energies E c _ m _ . The crosses denote the entrance channel potential with the Skyrme energy density 
funcitonal approach 5[ which is based on frozen density approximation. The open circle denotes the 
results of TDHF [32] at £7 c . m .=80MeV. The inserted sub- figures denote the the density distributions 
for this reaction at E c . m =80MeV and different relative distances, respectively. 

In this work, we write the nucleus-nucleus potential as a smooth function between V\ and 

V 2 , 

V b (R) = ^eric(s)V 2 + [1 - ^erfc( S )]F 1 (18) 
where erfc(s) is the complementary error function and 

with <5=lfm, Ai?=2fm. The obtained nucleus-nucleus potential in Eq.(18) approaches to 
V\ with the increase of the separation distance between two nuclei. On the contrary, Vb(R) 
approaches to V 2 with the formation of the di-nuclear system and the decrease of the distance 
between two nuclei. To study the dynamical nucleus-nucleus potential, we create 500 
reaction events for head-on collision of 40 Ca+ 40 Ca at several center-of-mass (cm.) energies 
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FIG. 4: (Color online) (a) Barrier height B for the reaction Ca+ Ca at different incident energies 
E c .m.- The horizontal dashed and solid lines indicate the barrier height of the entrance channel 
potential based on the frozen density approximation [5] and the lowest barrier Bq extracted from 
the fusion excitation function [8j], respectively, (b) The same as Fig. 1(b). 

ranging from E c _ m , = 52 to 140 MeV. For each event, we evolve the reaction system for 
a time of 700 fm/c. The distance between the projectile and target at the initial time is 
set to 30 fm for this reaction. The scattering events at t = 700 fm/c are not involved 
in the calculation of the nucleus-nucleus potential. Fig. 3 shows the obtained dynamical 
nucleus- nucleus potential at different incident energies E cm _. The corresponding density 
distributions at E c m =80MeV and different relative distances are also shown in the sub- 
figures. Fig. 4 shows the average fusion barrier height B for the reaction 40 Ca+ 40 Ca at 
different E cm . From the Fig. 3 and Fig. 4 we can see that the dynamical barriers depend 
on the incident energies. At energies around the Coulomb barrier, the dynamical barrier 
increases rapidly with the incident energy. With the further increase of the incident energy, 
the dynamical barrier approaches to the barrier height of the entrance channel potential 
(56.9 MeV) which is obtained with the Skyrme energy density functional together with the 



frozen density approximation [5]. This trend has also been found in [15l . |20j. When the 
incident energy decreases gradually and down to E c m ^ = 55 MeV, the height of dynamical 
barrier falls to 52.0 MeV. With lower incident energy, the height of the dynamical barrier 



10 



approaches to about 50 MeV which is close to the height of the lowest barrier Bq extracted 
from the fusion excitation function. In addition, it is encouraging that the obtained barrier 
height and the depth of the fusion pocket in this work are comparable with the results of the 
TDHF calculations |32J. The depth of the fusion pocket is about 25 MeV for this reaction 
system. 

At very short distances between two nuclei, it is thought that the Q value of the fusion 
system may provide information of the nucleus-nucleus potential. One commonly defines 
the excitation energy for a reaction from the expression 

— Qgg = -Ee.m. — E* , (20) 

where Q gg is the mass difference between the two initial nuclei and the combined system in 
its ground state. From Eq.(15) and Eq.(17), one can find that when the compound nucleus 
is well formed and the collective motion can be negligible one gets V(R = 0) = —Q gg , which 
is the result of the adiabatic nucleus-nucleus potential [l0|. For the reaction 40 Ca+ 40 Ca, we 
have — <5 99 =14.3 MeV. Actually, because the expression (20) is correct relative to the ground 
state of the composite system, it does not accurately describe the excitation energy relative 



lsion 



321 ]. the nucleus- nucleus 



32| at short distance does not 



to other intermediate transition states formed during the col . 
potential obtained from the ImQMD model and the TDHF 
exactly reach the value — Q gg since the composite system formed during the collision is far 
from a ground state. 

With the same approach we studied the nucleus-nucleus potential for the reactions 
48 Ca+ 208 Pb and 126 Sn+ 130 Te at energies above the Coulomb barrier. The corresponding 
values of — Q gg for these two reactions are 153.8 and 261.2 MeV, respectively. These two 
reactions lead to the same compound nucleus 256 No. Fig. 5 shows the calculated nucleus- 
nucleus potentials for these two reactions. The arrows denote the Bass barriers. From Fig. 
5 we can see that the obtained barrier heights are close to the corresponding Bass barriers. 
The depth of the fusion pocket (about 7 MeV) for 48 Ca+ 208 Pb becomes much shallower than 
that of 40 Ca+ 40 Ca (about 25 MeV) and the fusion pocket for 126 Sn+ 130 Te almost disappears, 
which indicates that quasi-fission could easily occur in heavy fusion process especially for 
the more symmetric systems. Furthermore, we find that the nucleus-nucleus potentials for 
the reactions 48 Ca+ 208 Pb and 126 Sn+ 130 Te at short distances are much higher than the value 
of — Q gg and even higher than the Coulomb barrier, which is quite different from the case 
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FIG. 5: (Color online) Dynamical nucleus-nucleus potentials for the reactions 48 Ca+ 208 Pb and 
126 Sn+ 130 Te at an incident energy £7 cm _=200MeV and £7 c . m ,=350MeV, respectively. The initial 
distance is set to 40 fm. The arrows denote the corresponding Bass barriers. 

of 40 Ca+ 40 Ca. These calculations indicate: 1) additional incident energy (so-called extra- 
push energy [2j) beyond the energy to overcome the Coulomb barrier may be required to 
form the compound nucleus for heavy fusion system and 2) the process of nucleon transfer 
between the projectile (like) and the target (like) could last for a period of time due to the 
appearance of the fusion pocket in di-nuclear system, which is the basic assumption of the 
di-nuclear system (DNS) model To see the fusion path, we also show the corresponding 
density distributions of the composite systems at t — 700/m in the inserted sub-figures. 
One can see that the corresponding compound nuclei are not well formed at t — 700/m 
for these two heavy fusion systems. The strongly deformed composite systems or called 
di-nuclear systems are formed at about t = 350 fm and can last hundreds even thousands 
fm/c for heavy fusion system, which is quite different from the case of light system such as 
40 Ca+ 40 Ca in which the spherical composite system is well formed at t — 700/m with the 
incident energies above the Coulomb barrier (see Fig.3). For 126 Sn+ 130 Te, the composite 
system tends to undergo quasi-fission or fission. In Fig. 6 we show the capture cross sec- 
tions of these two reactions. The solid and open circles in Fig. 6 (a) denote the experimental 
data of 48 Ca+ 208 Pb in Ref-jsiJ and 34], respectively. The solid curves denote the results 



of an empirical barrier distribution approach which is based on the Skyrme energy- density 
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FIG. 6: (Color online) Capture cross sections for the reactions Ca+ Pb and Sn+ Te. The 



33, 



34|]. The solid 



35(]. The 



open and solid circles in (a) denote the experimental data of 48 Ca+ 208 Pb 
curves denote the results of an empirical barrier distribution approach proposed in [5j, 
solid squares denote the results of the ImQMD model and the error bars denote the corresponding 
statistical error. 



functional together with the extended Thomas- Fermi (ETF) approximation p, |35|]. The 
solid squares denote the results of the ImQMD model with IQ2 and the error bars denote 
the corresponding statistical errors. For the reaction 48 Ca+ 208 Pb, the experimental data 
at energies above the Coulomb barrier can be reproduced acceptably well by the ImQMD 
model. Because the ImQMD model has difficulties to deal with the shell effects, the cap- 
ture cross sections of 48 Ca+ 208 Pb at sub-barrier energies can not be described well. For 
126 Sn+ 130 Te, the calculated capture cross sections with the ImQMD model are comparable 
with the results of the empirical barrier distribution approach. 

IV. SUMMARY 

In summary, the kinetic energies of a series of nuclei have been studied with the ImQMD 
model together with the extended Thomas-Fermi (ETF) approximation which gives accurate 
results for finite nuclear system, especially for the light and heavy nuclei. Furthermore, with 
the "ETF" approximation for the kinetic energies we have studied the dynamical Coulomb 
barrier of the reaction 40 Ca+ 40 Ca at different incident energies. The results show that the 
dynamical Coulomb barrier strongly depends on the incident energy. With the increase of 
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the incident energy, the dynamical Coulomb barrier increases gradually and approaches to 
the entrance channel potential which is based on the frozen density approximation. The 
height of dynamical Coulomb barrier decreases with the decrease of the incident energy and 
approaches to the lowest barrier extracted from the fusion excitation function. The behavior 
of nucleus-nucleus potential at short distances for heavy system is obviously different from 
that of light systems. For heavy fusion systems, the depth of the fusion pocket becomes 
much shallower and the nucleus-nucleus potential at short distances are higher than the 
adiabatic potential. The capture cross sections for 48 Ca+ 208 Pb and 126 Sn+ 130 Te have also 
been studied with the ImQMD model. The calculated results are comparable with the 
results of the empirical barrier distribution approach. A systematic study of heavy fusion 
systems, such as the calculation of the potential energy surface of the composite system in 
the fusion process as a function of mass-asymmetry and the distance between two nuclei is 
in progress with the shell effects being taken into account. 
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